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Observations support the view that feedback, in the form of radio outbursts from active nuclei in central galaxies, prevents 
catastrophic cooling of gas and rapid star formation in many groups and clusters of galaxies. Variations in jet power drive 
a succession of weak shocks that can heat regions close to the active galactic nuclei (AGN). On larger scales, shocks 
fade into sound waves. The Braginskii viscosity determines a well-defined sound damping rate in the weakly magnetized 
intracluster medium (ICM) that can provide sufficient heating on larger scales. It is argued that weak shocks and sound 
dissipation are the main means by which radio AGN heat the ICM, in which case, the power spectrum of AGN outbursts 
plays a central role in AGN feedback. 
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1 Introduction 

X-ray observations of cool core clusters and groups have 
shown that radio outbursts from central AGN have large 
impacts on their hot atmospheres (reviewed in McNamara 
& Nulsen 2007). Radio outbursts originating near the event 
horizons of supermassive black holes deposit energy on spa- 
tial scales eight orders of magnitude larger In the absence of 
a heat source, copious amounts of hot gas would be cooling 
to low temperatures and forming stars. Remarkably, pow- 
ers of AGN outbursts are comparable to the powers needed 
to stop the gas from cooling, signalling that the mechanical 
power output of the AGN is governed by feedback (Birzan 
et al. 2004; Dunn & Fabian 2006; Rafferty et al. 2006; Cro- 
ton et al. 2006; Sijacki & Springel 2006). Cooled or cooling 
gas can fuel AGN outbursts, while the outbursts heat the gas, 
affecting the fuel supply for subsequent outbursts. The high 
incidence of clusters with central cooling times less than 1 
Gyr (Hudson et al. 2010) also argues for radio mode AGN 
feedback. While many processes might heat the ICM and 
prevent the gas from cooling, without feedback, it is all but 
impossible to account for the many systems with very short 
cooling times, while essentially no clusters are undergoing 
catastrophic cooling (McNamara & Nulsen 2012). 

While the broad outline of the feedback cycle seems 
simple, very little of the detail is understood. In this article 
we focus on some processes that may heat the gas on scales 
comparable to the Bondi radius and larger Apart from the 
discussion of viscosity, much of this material has been re- 
viewed more thoroughly by McNamara & Nulsen (2012). 
We focus on the heating effects of weak shocks in section 2 
and sound waves in section 3. In section 4 we consider how 
these processes operate together. 
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2 Heating and outburst history 

2.1 Adiabatic uplift is ineffective 

To prevent the hot gas from cooling and forming stars, the 
key requirement is a heat source to replace heat lost by 
radiation. It should be emphasized that adiabatic uplift is 
ineffective at preventing gas from cooling. Extended fila- 
ments of cool, low entropy gas (e.g. Werner et al. 2010; 
Gitti et al. 201 1) and heavy elements shed by evolving stars 
in the central galaxy (e.g. Million et al. 2010; Kirkpatrick et 
al. 2011) make a strong case for gas uplift in the wakes of 
radio lobes. Lifting gas outward adiabatically, into regions 
where the atmospheric pressure is lower, will generally ex- 
tend its cooling time and so can help to delay the onset of 
cooling. However, for abundances and temperatures in the 
relevant range, the effect is modest. For example, for gas 
with solar abundances, starting at 3 keV and reducing the 
pressure adiabatically by a factor of ~ 88 would reduce its 
temperature to 0.5 keV, but only increase its cooling time by 
~ 36%. For gas with 0.5 solar abundances, the increase in 
cooling time is a little under a factor of 2, still well short of 
what is required to prevent the gas cooling in the long term 
(McNamara & Nulsen 2012). 

Furthermore, unless the uplifted gas mixes with it sur- 
roundings, raising its entropy, it is negatively buoyant and 
will fall back to where it came from in about one free-fall 
time. That is generally much shorter than the cooling time, 
so the effect of lifting the gas is transient. Metals shed by 
cluster central galaxies are more extended than their stars, 
showing that they have diffused outward from where they 
were shed (Rebusco et al. 2005). However, if all the uplifted 
gas mixed effectively the mean diffusion rate would be far 
too high, so most of the gas must fall back almost to where 
it originates. Much of the energy needed to lift the gas is 
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then converted to kinetic energy and dissipated in the gas, 
providing a channel for heating (Gitti et al. 201 1). 

2.2 Weak shocks 

The thermal energy of the gas within the volume V is 



^th-3/2 [ pdV, 
Jv 



(1) 



where p is its pressure. If the energy deposited by an AGN 
into this volume is comparable to or larger than Eth, then 
the fractional pressure increase in V must be large to ac- 
commodate the extra energy. That would cause the region 
affected by the jet to expand supersonically, driving shocks 
into its surroundings and rapidly extending the region af- 
fected by the outburst. Similarly, shocks will be formed if 
the jet power exceeds Ety^ divided by the sound crossing 
time of V. Thus, the region affected by an outburst must 
contain a thermal energy prior to the outburst that is sig- 
nificantly larger than the outburst energy and also larger 
than the jet power times the sound crossing time of the 
region. Otherwise, shocks are generated. Based on simula- 
tions, Morsony et al. (2010) have found that, under the in- 
fluence of cluster "weather" due to continuing infall, motion 
of substructures, etc., the radius of influence of the central 
AGN scales with the power of its jet, Pjot, as i?influcncc « 

1/3 

p. 'J. . The argument here suggests that this cannot be the 
whole story. 

A number of systems, such as MS0735.6H-7421 (McNa- 
mara et al. 2005) and NGC 5813 (Randall et al. 201 1), show 
symmetric, large-scale shocks fronts that do not appear to be 
affected significantly by the cluster weather. Furthermore, 
in NGC 5813 (Randall et al. 2011) and M87 (Forman et 
al. 2007), there is clear evidence of multiple cavities and 
shocks. The nested shocks seen in these systems almost cer- 
tainly require large and sustained variations in the power of 
the jet for their formation, consistent with other evidence of 
variations in jet power (e.g. Wise et al. 2007). Ripples in the 
Perseus ICM (Fabian et al. 2006) also indicate power varia- 
tions (though sound may be driven in other ways, e.g. Stem- 
berg & Soker 2009). Note that the small fluctuations seen 
at large distances from the cluster centre would have been 
considerably greater when they were launched in the region 
near NGC 1275. Power variations with shorter timescales 
and/or lower amplitude would be too weak to see now due 
to the combined effects of damping (section 3) and the de- 
crease in amplitude with increasing radius. It is notewor- 
thy that the two best observed systems with AGN outbursts, 
M87 and Perseus, show multiple weak shocks. Given the 
ubiquitous evidence of AGN variability on a broad range of 
timescales, this is probably the norm. As argued below, the 
power spectrum of AGN outbursts plays a critical role in 
AGN feedback by controlling the launching of shocks and 
sound waves. 

Although the heating effect of individual weak shocks is 
minor, their cumulative impact need not be. The changes in 
thermal and kinetic energy associated with a weak shock 



may be substantial, but mostly move on with the shock. 
A small entropy increases, AS", that is cubic in the shock 
strength (Landau & Lifshitz 1959), is all that remains. The 
heat equivalent of the entropy jump is AQ = T AS = 
E A\nK, where K = kT/nJ is the entropy index and 
A In X is the jump of In K in the shock. Expressed as a 
fraction of the gas thermal energy, AQ/E = AlnK. For 
the innermost shock in M87, at a radius of ~ 0.8 arcmin 
(3.7 kpc; Forman et al. 2007), the Mach number is ~ 1.38, 
giving an equivalent heat input of only AQ/E ~ 0.022. 
There is a second shock at about twice the radius and a 
third shock that is several times more energetic at a radius 
of ~ 3 arcmin. The shock spacings suggest that shocks of 
similar strength to the innermost shock are launched every 
~ 2.5 Myr, while the cooling time of the gas at 0.8 arcmin is 
~ 250 Myr. The ~ 100 shocks launched during one cooling 
time would add heat AQtot/E ~ 100 x 0.022 = 2.2, more 
than enough to replace the energy radiated. These numbers 
are indicative only, but they show that repeated weak shocks 
alone can prevent gas near the centre of M87 from cooling 
(Nulsen et al. 2007). A similar argument has been made for 
weak shock heating in NGC 5813 (Randafl et al. 201 1). 

Weak shocks can prevent gas cooling at the centres of 
two of the two nearest, best observed systems. The ripples in 
Perseus may well start as weak shocks launched from near 
the centre of NGC 1275, where they prevent the gas from 
cooling too. If the best observed systems are representative, 
weak shocks can be the primary channel for heating gas near 
the centres of all systems. Because shock strength declines 
with distance from the AGN and A In X depends steeply on 
shock strength, weak shocks become less effective at larger 
radii. However, as the shock strength decreases, sound dis- 
sipation increases in relative importance, probably taking 
over as the main heating channel, as outlined below. The 
efficiency of weak shock heating, measured by the fraction 
of shock energy converted to heat, is generally low, so that 
sound dissipation will often make a greater contribution to 
the total heating rate. 

3 Plasma viscosity and sound heating 

3.1 Plasma viscosity 

Although the magnetic pressure in the ICM is typically only 
^ 1% of the gas thermal pressure (Carilli & Taylor 2002), 
particle Larmor radii are ten or more orders of magnitude 
smaller than their mean free paths, keeping the particles tied 
rigidly to the magnetic field. As a consequence, transport 
processes in the ICM are poorly understood. This issue is 
worst for thermal conduction, since the heat flux depends 
on the structure of the magnetic field on scales comparable 
to the electron mean free path. The field may well vary on 
scales smaller than this, in which case the heat flux is not 
calculable in terms of local gas properties. Despite a great 
deal of work in this area, the role of thermal conduction in 
the ICM remains highly uncertain. 
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By contrast, viscous stresses are often determined lo- 
cally. Fluid motions readily push the dynamically insignif- 
icant field around and, being frozen in, the field generally 
varies as the fluid moves. In a collisionless plasma, the mag- 
netic moment, mv\/{2B), is conserved, where m is the 
particle mass, the strength of the magnetic field is B and 
the component of the particle velocity perpendicular to the 
field is v±. Thus changes in the magnetic field cause corre- 
sponding changes in the transverse kinetic energy, ■mv\/2, 
making the particle velocity distribution anisotropic. The 
proton-proton collision time is Tpp ~ 700(fcT)'^''^np ^ yr, 
where the temperature, kT, is in keV and the proton den- 
sity, rip, is in cm""^. Collisional relaxation can generally 
suppress most of the anisotropy caused by fluid motions, 
leaving only a small residual effect. Anisotropy in the ve- 
locity distribution implies anisotropy in the pressure, which 
can be expressed as the difference between the transverse 
and mean pressures, A = p^ — p. When collisional re- 
laxation is fast compared to the rate of change in B, i.e. 
Ti{\dB/dt\ <C B, where tj; is the ion-ion collision time, the 
pressure anisotropy is given by (e.g. Kunz et al. 2012) 



A = TiiPi ( bb : Vv - -V • V ) , 



(2) 



where pi is the ion pressure, b is the direction of the mag- 
netic field and v is the fluid velocity. The viscous stress ten- 
sor, which is minus the anisotropic part of the total stress 
tensor, is then 

T = A(3bb - 1), (3) 

where 1 is the 3x3 unit matrix. This is the Braginskii ( 1 965) 
form of the viscous stress tensor for a magnetized plasma. 
Anisotropy generated by motion parallel and perpendicular 
to the magnetic field simply reflects work done on or by 
the components of the particle velocities through the cor- 
responding particle pressures, p±^ or p|| = 3p — 2p^. In a 
uniform magnetic field parallel to the z axis, with gas mo- 
tions only parallel the field, T^z must have the same value 
as it does in the absence of a magnetic field. This requires 
TiiPi to equal the viscosity of an unmagnetized plasma, i.e. 
the "Spitzer" viscosity. Thus, although the viscous stresses 
depend on the direction of the magnetic field, they are simi- 
lar in magnitude to the stresses in the absence of a magnetic 
field. 

3.2 Sound heating 

Calculation of sound damping for a magnetized ICM is much 
the same as in the field free case. We assume that the unper- 
turbed magnetic field is uniform and parallel to the z axis. 
For sound with wavevector k — {kx,ky,kz), the damping 
rate of the amplitude is 



j//c^ 1 



m 



(4) 



where i^ = Tapi/p is the kinematic viscosity and p is the 
gas density. For sound waves travelling parallel to the mag- 
netic field (k^ — k^), this is identical to the damping rate 



with no magnetic field, 2z^fc^/3, as we should expect. For 
sound travelling perpendicular to the field (kz — 0), the 
damping rate is vk^ /&, a factor of 4 smaller. Note that for 
k1 ~ k"^ I'i, i.e. for sound waves inclined at ~ 54.7° to the 
magnetic field, the damping rate is zero. In this direction 
sound waves have the same effect on p±^ and p||, so they 
generate no anisotropy, hence no viscous damping. 

Due to cluster weather, the field is likely to be fairly 
chaotic, with structure on scales down to 10 kpc or less 
(e.g. Govoni et al. 2010). If the field is isotropic on aver- 
age, we can average the damping rate (4) over the sphere to 
obtain the mean value of 2i^fc^/15, one fifth of the damping 
rate with no magnetic field. Fabian et al. (2005) have shown 
that viscous damping of sound is a viable mechanism to pre- 
vent gas cooling in the Perseus cluster and Abell 2199. Their 
models used a viscous damping rate that is one half of the 
mean rate calculated here. As they discuss, the power spec- 
trum of the AGN outbursts is critical for the heating rate, 
since it controls the spectrum of sound generated by out- 
bursts and the dissipation rate is sensitive to the frequency. 
Corresponding to the average damping rate, the dissipation 
length for sound power is 

~' ^hL__\ ( ^^ V' 

.0.03cm--V UkcVJ 



120 



kpc, (5) 



where the angular frequency unit, Wg — 2.36 x 10^^'* s^^, 
gives a wavelength 10 kpc for a gas temperature of 5 keV 
(cf. Fabian et al. 2005). This is comparable to the size of 
cluster cool cores. If the thermal conductivity is close to 
its field free value, conduction would boost the dissipation 
rate significantly. However, viscous dissipation alone ap- 
pears sufficient to heat the ICM in Perseus and Abell 2199. 
Schekochihin et al. (2005) note that, when the anisotropy 
falls outside the range — 2pb/3 < A < pb/3, for ps = 
B'^ /{8tt), the plasma is prone to firehose or mirror insta- 
bilities that grow much faster than viscous damping. For 
Pb ^ P, as in clusters, this condition is more restrictive 
than the one preceding equation (2). It requires ljthA < 
0.5pb/p, where A is the amplitude of the fractional den- 
sity fluctuations. If the anisotropy did get too large, sound 
damping could be much faster than the viscous rate and 
sound heating even more effective — probably too effec- 
tive, although some sound energy would then be converted 
to magnetic field rather than being thermalized (Schekochi- 
hin et al. 2005). 

4 Discussion 

Variations in the power of the jet launched from the AGN 
initiate a succession of weak shocks that originate from some- 
where near the Bondi radius (disregarding the small amount 
of gas within the Bondi sphere). The distribution of gas 
close to the radio lobes is generally far from spherical and 
the shock strength varies over the surface of spheres (e.g. For- 
man et al. 2005; Mendygral et al. 2012). However, shock 
strength generally decreases radially. As the shocks weaken. 
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at some point they can be regarded as sound waves. Al- 
though the power spectrum of the sound is far from monochro- 
matic, for the purpose of discussion here, we assume that 
there is a representative frequency, to, where 2tt/uj is typi- 
cal of the time interval between shocks. 

Heating by weak shocks (section 2.2) depends on the 
fractional pressure jump, Sp/p, as {5p/pY, whereas sound 
damping heats the gas at a rate proportional to the square 
of the sound amplitude. As a result, shock heating is more 
significant at small radii where the pressure disturbance is 
larger, but as the disturbance decreases with increasing ra- 
dius, sound dissipation will overtake it at some point. If the 
mean kinematic viscosity is z^/5 (section 3.2) and the frac- 
tional amplitude of the sound pressure disturbance is related 
to the pressure jump by Ap = 0.56p/p, the ratio of the 
shock heating rate to the sound heating rate is (McNamara 
& Nulsen 2007) 
2s2 /5p^ 



P 



12.6 



kT 



-3/2 



-] (6) 

in the units of equation (5). If the entropy or the effective 
frequency is high (for wavelengths comparable to the pro- 
ton mean free path or smaller), sound damping might ri- 
val shock heating, even for Sp/p ^ 1, but low entropies in 
cool cores favour weak shock heating. Weak shock heating 
probably dominates close to the Bondi radius, while sound 
damping is more effective on larger scales. For both pro- 
cesses, the power spectrum of AGN outbursts plays a cen- 
tral role in determining the heating rate. Clearly, we need to 
understand how that is governed. 

5 Conclusions 

Low rates of gas cooling and star formation in the central 
galaxies of many groups and clusters are best explained by 
feedback from radio AGN. This view is supported by obser- 
vations of the impacts of radio outbursts in hot atmospheres. 
Jet powers vary significantly on a wide range of timescales, 
causing a succession of shocks to be launched from near the 
radio lobes. The shocks weaken into sound waves on larger 
scales. For the best observed, nearby systems, weak shocks 
provide sufficient heat near the AGN to prevent gas from 
cooling. 

In contrast to thermal conduction, the Braginskii viscos- 
ity generally provides a well-defined, local, viscous stress 
that only depends on the direction of the magnetic field in 
the weakly magnetized ICM. It determines a viscous damp- 
ing rate for sound that is a factor of five smaller than that for 
unmagnetized plasma when the field is isotropic on aver- 
age. This level of sound damping is sufficient to prevent gas 
from cooling throughout the cool core of the Perseus cluster. 
The combination of weak shock heating on small scales and 
sound damping on larger scales plausibly provides the pri- 
mary means by which AGN energy heats the ICM and pre- 
vents gas from cooling. If so, the power spectrum of AGN 



outbursts plays a central role in AGN feedback and needs to 
be better understood. 
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